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Role of the isospin asymmetry in nuclei and neutron stars, with an 
emphasis on the density dependence of the nuclear symmetry energy, is 
discussed. The symmetry energy is obtained using the isoscalar as well 
as isovector components of the density dependent M3Y effective interac- 
tion. The constants of density dependence of the effective interaction are 
obtained by reproducing the saturation energy per nucleon and the satu- 
ration density of spin and isospin symmetric cold infinite nuclear matter. 
Implications for the density dependence of the symmetry energy in case of 
a neutron star are discussed, and also possible constraints on the density 
dependence obtained from finite nuclei are compared. 

PACS numbers: 21.65.+f, 26.60,+c, 97.60.Jd, 21.30.Fe 

1. Introduction 

The symmetry energy is an important quantity in the equation of state 
of isospin asymmetric nuclear matter. This currently unknown quantity 
plays a key role to the understanding of the structure of systems as diverse 
as the neutron rich nuclei and neutron stars. Theoretical studies based on 
microscopic and many-body calculations and phenomenological approaches 
predict various forms of the density dependence of the symmetry energy. In 
general, they can be broadly classified into two different forms. One, where 
the symmetry energy increases monotonically with increasing density ('stiff' 
dependence) [lj and the other, where the symmetry energy increases [2] ini- 
tially up to normal nuclear density or slightly beyond and then decreases 
at higher densities ('soft dependence'). Determination of the exact form of 
the density dependence of the symmetry energy is of utmost importance 
for studying the structure of neutron rich nuclei, and studies relevant to 
astrophysical problems, such as the structure of neutron stars and the dy- 
namics of supernova collapse [3]. For example, a 'stiff' density dependence 
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of the symmetry energy is predicted to lead a large neutron skin thickness 
compared to a 'soft' dependence and can result in rapid cooling of a neutron 
star, and larger neutron star radius, compared to a soft density dependence. 

In view of rather large differences between the various calculations of 
the symmetry energy present even at the subsaturation densities, the ques- 
tion arises whether one can obtain empirical constraints from finite nuclei. 
Since the degree of isospin diffusion in heavy-ion collisions at intermediate 
energies is affected by the stiffness of the nuclear symmetry energy, these re- 
actions, therefore, can also provide constraints on the low energy behaviour 
of the nuclear symmetry energy pQ. However, the high density behaviour re- 
mains largely undetermined since the masses and corresponding radii of the 
neutron stars are not measured whereas they can be obtained theoretically 
by solving Tolman-Oppenheimer-Volkov equation. At present there exist 
only indirect indications such as the neutron star cooling process. Recently 
search for the experimental signatures of the moderately high density be- 
haviour of the nuclear symmetry energy has been proposed [2] theoretically 
using several sensitive probes suh as the 7r~ to tt + ratio, tranverse collective 
flow and its excitation function as well as the neutron-proton differential 
flow. In the present work, the nuclear symmetry energy is calculated theo- 
retically using the isoscalar and isovector components of M3Y-Reid-Elliott 
effective interaction supplemented by a zero range pseudo-potential [4] along 
with the density dependence and its high density behaviour is explored. 



2. The nuclear symmetry energy 

The nuclear symmetry energy (NSE) E sym (p) is defined as the energy re- 
quired per nucleon to change the spin and isospin symmetric nuclear matter 
(SNM) to the spin symmetric pure neutron matter (PNM). In the present 
work the NSE is obtained using the density dependent M3Y interaction 
supplemented by a zero-range pseudo-potential (DDM3Y) as the effective 
interaction between two interacting nucleons inside nuclear medium. The 
M3Y interaction was derived by fitting its matrix elements in an oscillator 
basis to those elements of the G-matrix obtained with the Reid-Elliott soft- 
core NN interaction. The ranges of the M3Y forces were chosen to ensure a 
long-range tail of the one-pion exchange potential as well as a short range 
repulsive part simulating the exchange of heavier mesons [S] . The real part 
of the nuclear interaction potential obtained by folding in the density distri- 
bution functions of two interacting nuclei with the DDM3Y effective inter- 
action was shown to provide good descriptions for medium and high energy 
a and heavy ion elastic scatterings [U [6] . The zero-range pseudo-potential 
represented the single-nucleon exchange term while the density dependence 
accounted for the higher order exchange effects and the Pauli blocking ef- 
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fects. The real part of the proton-nucleus interaction potential obtained by 
folding in the density distribution function of interacting nucleus with the 
DDM3Y effective interaction is found to provide good descriptions of elastic 
and inelastic scatterings of high energy protons [7] and proton radioactivity 

The central part of the effective interaction between two nucleons 1 and 
2 can be written as [I] 



v 12 (s) = v 00 (s) +voi(s)T 1 .T2 + v w (s)a 1 .a 2 +vn{s)ai.a2 n.r 2 , (1) 

where t\,t 2 are the isospins and ai,cr 2 are the spins of nucleons 1,2 and 
s is the distance between them. In case of SNM, only the first term, the 
isoscalar term, contributes whereas for the isospin asymmetric-spin symmet- 
ric nuclear matter only first two terms, the isoscalar and the isovector terms, 
contribute and for the spin-isospin asymmetric nuclear matter all the four 
terms of Eq.(l) contribute. Considering only the isospin asymmetric-spin 
symmetric nuclear matter, the neutron-neutron, proton-proton, neutron- 
proton and proton-neutron interactions, viz. v nn ,v P p,v np and v pn respec- 
tively, can be given by the following : 

Vnn = V pp = VQO + V 1, V np = V pn = V ~ «01- ( 2 ) 

The general expression for the density dependent effective NN interaction 
potential is written as [8] 



v o(s,p,e) = t^ Y (s,e)g(p,e), v 0l (s, p, e) = t™ iY (s, e)g(p, e), (3) 

where e = E/A is the energy per nucleon and the isoscalar ioo 3 ^ an d the 
isovector iof 3y components of M3Y interaction potentials [4J supplemented 
by zero range potentials are given by the following : 



t^(s, e) = 7999^=M - 2134^g^ - 276(1 - ae)5(s) (4) 

and 



4{ 3¥ (s, e) = -4886^^1 + 1176^1^ + 228(1 - ae)5(s), (5) 

4s 2.5s 

where the energy dependence parameter a=0.005/MeV. The zero-range po- 
tentials of Eqs.(4,5) represent the single-nucleon exchange term. The density 



4 



NSE printed on February 9, 2008 



dependent part appearing in Eq.(3) [9] has been taken to be of a general 
form 

g(p,e)=C(l-(3(e)p n ) (6) 

which takes care of the higher order exchange effects and the Pauli blocking 
effects. This density dependence changes sign at high densities which is of 
crucial importance in fulfilling the saturation condition as well as giving 
different incompressibility Kq values for SNM in ground state with different 
values of n for the nuclear equation of state (EOS) [9j. The value of the 
parameter n = 2/3 was originally taken by Myers in the single folding 
calculation [10]. The isospin asymmetry X can be conveniently defined as 

X = ^E, P = Pn + P P , (7) 

Pn + Pp 

where p n , p p and p are the neutron, proton and nucleonic densities respec- 
tively. The asymmetry parameter X can have values between -1 to +1, 
corresponding to pure proton matter and pure neutron matter respectively, 
while for SNM it becomes zero. For a single neutron interacting with rest of 
nuclear matter with isospin asymmetry X, the interaction energy per unit 
volume at s is given by the following : 



PnVnn{s) + PpV n p(s) = p n [v o(s) + V 01 (s)] + p p [v o(s) - V 01 (s)] 

= [v 00 (s) + v 01 (s)X]p, (8) 

while in case of a single proton interacting with rest of nuclear matter with 
isospin asymmetry X, the interaction energy per unit volume at s is given 
by the following : 



PnVpn(s) + PpVpp(s) = p n [voo(s) - V l{s)] + p p [v 00 (s) + V 01 (s)\ 

= [voa(s) - v 01 (s)X]p. (9) 

Summing the contributions for protons and neutrons and integrating over 
the entire volume of the infinite nuclear matter and multiplying by the 
factor i to ignore the double counting in the process, the potential energy 
per nucleon e pot can be obtained by dividing the total potential energy by 
the total number of nucleons, 

_ g(p,e)pJ v 

Cpot — ^ , UUJ 

where 
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Jv = Jvoo + X 2 J v01 = J J J K 3Y (s, <0 + t& 3Y (s, e)X 2 ]d 3 s. (11) 

Assuming interacting Fermi gas of neutrons and protons, the kinetic 
energy per nucleon t^ in turns out to be 



.^[flW, f W - , 01^(1^! ], (12) 

where m is the nucleonic mass equal to 938.91897 MeK/c 2 and &f, which 
becomes equal to Fermi momentum in case of the SNM, is given by the 
following : 

k% = 1.5vrV (13) 

The two parameters of Eq.(6), C and /3, are determined by reproduc- 
ing the saturation conditions. It is worthwhile to mention here that due 
to attractive character of the M3Y forces the saturation condition for cold 
nuclear matter is not fulfilled. However, the realistic description of nu- 
clear matter properties can be obtained with this density dependent M3Y 
effective interaction. Therefore, the density dependence parameters have 
been obtained by reproducing the saturation energy per nucleon and the 
saturation nucleonic density of the cold SNM. 

The energy per nucleon e = ekin + tpot obtained for the cold SNM for 
which X = is given by the following : 

= 3h 2 k 2 F g(p, e)pJ v00 
6 [ 10m 1 " 2 ' ' [ ' 

where J v oo(e) = J J J tQQ 3Y (s,e)d 3 s represents the volume integral of the 
isoscalar part of the M3Y interaction supplemented by the zero-range po- 
tential. The Eq.(14) can be rewritten with the help of Eq.(6) as 

t 3h 2 k? F , , [ P J v o C(l-Pp n ) i 

e{p) = [ To^T ] + [ 2 ] (15) 

and differentiated with respect to p to yield equation 

The equilibrium density of the cold SNM is determined from the saturation 
condition = 0. Then Eq.(15) and Eq.(16) with the saturation condition 
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can be solved simultaneously for fixed values of the saturation energy per 
nucleon eo and the saturation density po of the cold SNM to obtain the 
values of the density dependence parameters (5 and C. Density dependence 
parameters f3 and C, thus obtained, can be given by the following : 

= (17) 
P [(3n + 1) - (n + l)p] ' 1 } 

where 
and 

k Fo = [1.5tt 2 po] 1/3 , (19) 

Q — 1 fQ 1 (20) 

5mJ,x)oPo[l - (n + 1)Ppq] ' 

respectively. It is quite obvious that the density dependence parameter (3 
obtained by this method depends only on the saturation energy per nucleon 
eo, the saturation density po and the index n of the density dependent part 
but not on the parameters of the M3Y interaction while the other density 
dependence parameter C depends on the parameters of the M3Y interaction 
also through the volume integral J v qq. 

The incompressibility Kq of the cold SNM which is defined as 

can be theoretically obtained using Eq.(13), Eq.(16) and Eq.(21) as 

5m 2 

Since the product J v qqC appears in the above equation, a cursory glance 
reveals that the incompressibility Kq depends only upon the saturation 
energy per nucleon eo, the saturation density po and the index n of the 
density dependent part of the interaction but not on the parameters of the 
M3Y interaction. 

The energy per nucleon for nuclear matter with isospin asymmetry X 
can be rewritten as 
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3h 2 k 2 F 



}F(X) + ( 




10m 



3h 2 k 2 F 



]F(X)-[f][ 



p j v h 2 k 2 Fo (i-(3 P n ) 



(23) 



10m 



po Jvoo 5m[l - (n + l)f3p%] 



where J v = J v qq + X 2 J v q\ and J v oi(e) = J J J iof 3y (s, e)<i 3 s represents the 
volume integral of the isovector part of the M3Y interaction supplemented 
by the zero-range potential. 

3. Density dependence of nuclear symmetry energy and 



The nuclear symmetry energy (NSE) E sym (p) which is the energy re- 
quired per nucleon to change the SNM to pure neutron matter (PNM) is 
defined as 



The calculations have been performed using the values of the satura- 
tion density po = 0.1533/m -3 [11] and the saturation energy per nucleon 
eo = — 15.26MeU |12j for the SNM obtained from the co-efficient of the vol- 
ume term of Bethe-Weizsacker mass formula which is evaluated by fitting 
the recent experimental and estimated atomic mass excesses from Audi- 
Wapstra-Thibault atomic mass table [13J by minimizing the mean square 
deviation. For a fixed value of /?, the parameters a and C can have any pos- 
sible simultaneous values as determined from SNM. Using the usual values of 
a = 0.005MeU _1 for the parameter of energy dependence of the zero range 
potential and n = 2/3 [9], the values obtained for the density dependence 
parameters C and f3 and the SNM incompressibility Kq are 2.07, 1.668/m 2 
and 293.4 MeV respectively. The general theoretical observation by Colo' et 
al. is that the non-relativistic [2] and the relativistic [15] mean field models 
predict for the Kq values which are significantly different from one another, 
namely « 220-235 MeV and rs 250-270 MeV respectively. Considering the 
uncertainties in the extractions of eo [12] and po values from the experi- 
mental masses and electron scattering, present non-relativistic mean field 
model estimate for Kq is rather close to the relativistic mean field model 
estimates. However if one uses n = 1/3 instead of n = 2/3 then the value 
of Kq comes out to be 226.1 MeV which is close to the lower limit of the 
other non-relativistic calculations. Using Eq.(15) and Eq.(23), the NSE is 
given by 



constraints from finite nuclei 



Esyrn(p) = e(p, 1) - e(^,0). 



(24) 
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Fig. 1. The energy per nucleon e = E/A of SNM (spin and isospin symmetric nuclear 
matter), PNM (pure neutron matter), NSE (nuclear symmetry energy E sym ) are 
plotted as functions of p/po for the present calculations using DDM3Y interaction 
and the saturation energy per nucleon equal to -15.26 MeV. 



E sym {p) = (2 2 / 3 - l)^E° F (-) 2/3 + - /W.oi, (25) 

£> Po 2 

where the Fermi energy E% = f or the SNM at ground state. The 

first term of the right hand side is the kinetic energy contribution with 
density dependence of /? 2//3 whereas the second term arising due to nuclear 
interaction has a density dependence of the form of a%p + d2p n+1 with a\ 
and ci2 as constants with respect to density and values of n are limited to 
1/3 to 2/3 for any reasonable values of incompressibility. It is interesting 
to note that the effective interaction which provides unified description for 
radioactivity, nuclear matter and nuclear scattering with n = 2/3 predicts 
the density dependence for the NSE to be of the form aop 2 / 3 + ciip + d2p 5 / 3 - 
In Fig.-l plots of E/A for SNM, PNM and NSE as functions of ^ are 
shown for n = 2/3. The density dependence of the NSE at subnormal den- 
sity from isospin diffusion [I] in heavy-ion collisions at intermediate energies 
has an approximate form of 31. 6[-^-] MeV. In Fig. -2 this low energy be- 
haviour of NSE « 31.6[^] L05 MeV is plotted along with NSE obtained 
using Eq.(25) for two extreme values of n = 1/3 and 2/3. At a density 
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Fig. 2. The NSE (nuclear symmetry energy E sym ) as a function of p/po for the 
present calculations using DDM3Y interaction is compared with the low energy 
behaviour extracted from isospin diffusion in heavy-ion collisions. 



around /> ~ po both curves gives the same value of the NSE and at sub- 
normal densities both forms are very close. At higher energies the present 
NSE using DDM3Y interaction peaks at p ~ 1.8po an d becomes negative 
at p ~ 4po- A negative NSE at high densities implies that the pure neutron 
matter becomes the most stable state. Consequently, pure neutron matter 
exists near the core of the neutron stars and since E/A for PNM is always 
positive, it is unbound by the nuclear force but bound due to gravitational 
attraction. 



4. Cooling of neutron stars and the nuclear symmetry energy 

The knowledge of the density dependence of nuclear symmetry energy is 
important for understanding not only the stucture of radioactive nuclei but 
also many important issues in nuclear astrophysics, such as nucleosynthesis 
during presupernova evolution of massive stars and the cooling of protoneu- 
tron stars. A neutron star without neutrino trappings can be considered 
as a n,p,e matter consisting of neutrons (n), protons (p) and electrons (e). 
The neutrinos do not accumulate in neutron stars and the (3 equilibrium 
proton fraction xp [= p p /p] is determined by [16] 



ftc(3vrV/3) 1/3 = ±E sym (p)(l - 2xp). 



(26) 
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Fig. 3. The f3 equilibrium proton fraction calculated with NSE (nuclear symmetry 
energy) obtained using DDM3Y interaction is plotted as a function of p/po- 



The (3 equilibrium proton fraction is therefore entirely determined by the 
NSE. The (3 equilibrium proton fraction calculated using the present NSE 
is plotted as function of in Fig. -3. The maximum of xp ~ 0.048 occurs 
at p « 1.4/Oo and goes to zero at p ps 3.9po for n = 2/3. The NSE extracted 
from the isospin diffusion in the intermediate energy heavy- ion collisions, 
having the approximate form of 31. 6[-^] 05 MeV, provides a monotonically 
increasing (3 equilibrium proton fraction and therefore can not be extended 
beyond normal nuclear matter densities. Present calculation, using NSE 
given by Eq.(25), of the (3 equilibrium proton fraction forbids the direct 
URCA process since the equilibrium proton fraction is always less than 1/9 
|16j which is consistent with the fact that there are no strong indications that 
fast cooling occurs. Although SNM incompressibilty is slightly on the higher 
side, yet the present calculations provide a rather 'soft' nuclear symmetry 
energy. Using n = 1/3 for which SNM incompressibility is 226.1 MeV and 
close to the lowest acceptable limit, results are very much similar for the (3 
equilibrium proton fraction excepting that it goes to zero at p 5.1po- 



5. Summary and conclusion 



In summary, we have calculated the symmetry energy using the isoscalar 
and the isovector components of M3Y effective NN interaction. The low 
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density behaviour of the symmetry energy is found to be consistent with the 
nuclear symmetry energy extracted from the isospin diffusion in heavy-ion 
collisions at intermediate energies. Although the nuclear incompressibility 
is on the higher side, yet the present calculations provide a 'soft' nuclear 
symmetry energy. The calculated [3 equilibrium proton fraction forbids the 
direct URCA process which is consistent with the fact that there are no 
strong indications that fast cooling occurs. Present theoretical calculations 
provide density dependence of the nuclear symmetry energy that can be 
verified by the experimental signatures in several sensitive probes suh as 
the 7r~ to 7T+ ratio, tranverse collective flow and its excitation function as 
well as the neutron-proton differential flow. 
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